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ABSTRACT: The s-wave tt — tt scattering lengths a^(T) at finite temperature T 
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^ ■ and isospin / = 0, 2 are calculated within the SU(2) Nambu-Jona-Lasinio model. o?{T) 

O , displays a singularity at the Mott temperature Tm, defined as m.^{TM) = '2mq{TM), while 

^ i a^{T) is singular in addition at the lower temperature T^, where mcr{Td) = 2m.j^{Td), 

0\ [ rricr and rriT^ being the masses of the a and tt mesons, respectively. Numerically we find 

^ ' Td = 198MeV and Tm = 215MeV. We speculate on possible experimental consequences. 
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The scattering of pion by pions, tt+tt — > tt+tt is one of the most fundamental hadronic 
processes of QCD at a mesonic level. This is so for several reasons. In particular, this 
scattering involves only the lightest pseudoscalar modes of the theory, that occupy a special 
role as being the Goldstone particles associated with the almost exact SUl{2) x SUr{2) 
symmetry that is observed in the particle spectrum. As such, tt — tt reactions provide a 
direct link between the theoretical formalism of chiral symmetry and experiment. This 
is exemplified in the many existing studies of tt — tt scattering using chiral Lagrangians 
[1-7] that endeavour to calculate scattering amplitudes and scattering lengths. References 
[1-3] introduce [1] and use [2,3] chiral perturbation theory (ChPT), while [4-6] calculate 
scattering lengths in the Nambu-Jona-Lasinio (NJL) model [8,9]. Both s and p wave 
scattering lengths are well reproduced in these models. Essentially, the Weinberg limit 
[10] can be recovered in all cases, and corrections to it, due to processes such as pion 
rescattering are accounted for, for example in [2]. 

To date, investigations performed have all been at zero temperature. Our interest 
here, however, is to examine tt — tt scattering at finite temperature, but at zero baryon 
density, with a view to describing the baryon free central rapidity region occurring in ultra 
relativistic heavy ion collisions, in which pions are copiously formed. In this calculation, 
we focus in particular on those temperatures in the vicinity of the phase transition that 
takes the system from a chirally symmetric state at high temperature to a phase in which 
chiral symmetry is broken. In our study, we will restrict ourselves to the SU(2) version of 
the NJL model, 

C{x) = ij{x){i jd - mo)ip{x) + G[{iP{x)il;{x)f + {ij{x)n5rip{x))\ (1) 

where G is a coupling strength of dimension [Mass]~^, and mo is the current quark 
mass. Once again, we recall that although this model lacks confinement and is non- 
renormalizable, it gives a good and transparent description of the low energy meson sec- 
tor. One simply regulates divergent integrals with a cutoff. In our case, we will do so 
by restricting the 3-momentum of a quark or antiquark, \p\ < A. At a finite temperature 
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Tc, this model displays a second order phase transition in the chiral limit (mo = 0) from 
a broken to a chirally symmetric state. At finite values of the current quark mass mo, 
the phase transition is washed out [11]. The temperature range around Tc, i.e. 150 — 200 
MeV, is of particular importance for two reasons: (i) current experiments at CERN and 
AGS energies appear to produce a hadron plasma in this temperature range, and (ii) in 
the neighbourhood of second order phase transitions, divergences in measurable quantities 
often occur, such as the phenomenon of critical opalescence [12]. In our case, we might 
expect that certain cross sections display a rapid temperature dependence. It is therefore 
of interest to us to study the properties of tt — tt scattering in the vicinity of this phase 
transition. Our calculation of tt — tt scattering is performed by keeping only the diagrams 
from the lowest order expansion in l/A^c, with Nc the number of colors [13]. 

To establish our notation, we briefly summarize the basic formalism of pion-pion 
scattering. The reader may refer to Refs.[14]-[16] for a more complete discussion. Since 
three isospin channels are available for the process tttt — > tttt, the invariant scattering 
amplitude can be written in terms of three unknown functions, 

< CPc, dpd\T\apa:, bpb >= %b;cd = ^(s, t, u)5abScd + B{s, t, u)5acSbd + C{s, t, u)5adSbc- (2) 

Here a,b and c,d are the isospin labels of the initial and final states respectively, and s,t 
and u are the usual Mandelstam variables, s = {Pa+Pb)^^ t = {Pa—PcY' and u = {Pa—PdY- 
Perfect crossing symmetry enables one to express Tab-cd in terms of one amplitude alone, 
since under the exchange s <-> t, A{s,t,u) = B{t,s,u), and under s <-> tt, A{u,t,s) = 
C{s,t,u). It is a standard exercise to project out amplitudes of definite total isospin 
/, which we denote as {s^t^u). In the limit of scattering at threshold, \/i = 2m7r, 
t = tt = 0, the approach the scattering lengths a^, and one finds the relation 

To lowest order in 1/A^c) the invariant amplitude Tab-cd is calculated from the box and cr- 
propagation diagrams shown in Fig.l. Since these need only be evaluated at threshold, we 
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have pI = = = s/4. This choice of the pion momentum does not restrict generahty 
at T = since the system is Lorentz invariant. For T 7^ 0, the heat bath defines a special 
system and the specification pj = = s/A restricts our results to the scattering of pions 
whose cm. system is at rest in the heat bath. 

Box diagram: The box diagram of Fig.l is generic for all diagrams which can be 
constructed with internal quark and antiquark lines. From crossing relations, there are 
three such possibilities, which we label 1,2 or 3, and which correspond to the diagram 
shown, and those generated by s ^ t and s ^ u exchanges, respectively. Here, we follow 
the notation and calculation of Ref.[6]. A direct evaluation of the diagrams leads one to 
the result 

{ri)ah;cd = {SabScd + SacSbd - 5ad5bc)[^N^N fig^J[I (0) + I{p) - p^K{p)] 

{T2)ab;cd = {^abScd ' ^aJbd ' ^adhc) [^N.Nftg^^J [/(O) + lip) - p^Kip)] (4) 

{r3)ab;cd = i-SabScd + SacSbd + SadSbc)[SNcN fig^jl (0) + p'' L{p) /2 - 2p'K{p)] 

in terms of the integrals /(p), K{j>) and L{j>), 

m = / 
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^^^^ / (27r)4 -m2]2[(A;+p)2 -m2] 
L{p) = 



(27r)4 - to2]2[(A; + p)2 _ ^2)2 " 

For T = 0, the integrals in Eq.(5) have been evaluated by [6] to obtain the scattering 
lengths. At finite temperatures, however, one has J (i^/c/(27r)'^ — > ^ J2n I d'^k/{27v)^. Here 
P is the inverse temperature, and the sum on n runs over the Matsubara fermion frequen- 
cies i(jJni = (2n + l)7r/P, n = 0,±1,±2- • •, that occur in the 4-vector k = (iounik). 
Simple analytic expressions can be obtained for these functions, and they are listed in the 
appendix. 

a -propagation diagram: The diagram with intermediate a propagation shown in Fig. 
1 can be simply expressed in terms of its component parts. Again two further diagrams 
are generated by s <-> t and s <-> tt exchange. Defining the a — tt — tt vertex as 

r---(y,j9) = N,Nf J ^Tr[S{k+p)^,S{k)^,S{k+p')], (6) 



one may easily verify that the required choices of momenta for s and t channel graphs lead 
to 

r^'^'^(p, -p) = -8NcNfmI{p) (s - channel) 

and (7) 

r^'^'^(p, p) = -8NcNfm[I{0) - p^Kip)] (t - channel). 
Required also are the meson scattering amplitudes and D^, defined as 

-'^"f*) = i-2GnM(^ - 

for M = (T,7r. Since it can be shown that l-2GIlM{k) = mQ/m + AiGNcNf{k'^ -e\^)I{k), 
with = 0, €0- = 4m^, and also that = —mQ/mAiGN(.NfI{m.j^), one may express the 
functions in (8) in terms of I{p) only. One has 

[-iDM{k)]-^ = 2N,Nf[{k' - el,)I{k) - m^/(m,)]. (9) 

Now one may construct the amplitudes for the scattering diagrams. One has 

{rA)ab-cd = 5a65ed^^,,[r---(p, -p)]^D,{2p) (10) 

in the s-channel, and s ^ t and s u-exchanges give rise to 

{T5)ab;cd = (5„Ad + 5ad<^6c)4,,[r---(p,p)]2D,(0). (11) 

Finally, expressions (4), (10) and (11) require the temperature dependent pion-quark 
coupling QTrgq, and which, using the same arguments that led to Eq.(9), is given as 

9^,\ = -N'm + m - mlK{p)]\ 

It is now a simple matter to extract the functions A{s,t,u), B{s,t,u) and C{s,t,u) 
by summing the amplitudes {%)ab;cd of (4), (10) and (11), and identifying A,B, and C as 
the coefficients of SabScd, ^ac^bd and Sad^bc respectively. Prom these, one can construct the 
r^, and via Eq.(3), the s-wave scattering amplitudes can be obtained. One finds 

r° = GTi - Ta + 3% + 2% 

= (12) 
= 2% + 2% 



where the % are the functions in Eqs.(4), (10) and (11), here to be understood to be 
denuded of their isospin cofactors. One observes that is identically zero, as it must be 
for s- waves. 

Individual contributions to the scattering amplitude are shown as a function of tem- 
perature in Fig. 2. For this calculation, a standard set of parameters has been used, viz. 
A = 631MeV, mo = 5.5MeV and G = 5.514GeV~^ that lead to a value of the current 
quark mass m = 339MeV, m-,^ = 139MeV and /tt = 93.3MeV. For mo = 0, these values 
of A and G lead to a phase transition at = 195MeV. One observes that the T- matrix 
amplitudes are constant at low temperatures, and are approximately equal in magnitude 
at T = 0. Structure appears in the amplitudes 7^ at high values of the temperature. Ti and 
Ta are negative and divergent at the Mott temperature Tm, defined as the temperature 
at which the pion mode enters into the continuum, thereby dissociating into constituent 
quark and antiquark, i.e. 

m^(TM) = 2mq{TM). (13) 

The amplitude is positive and divergent at this temperature. On the other hand 74 
diverges at the temperature T^, at which the cr-meson can dissociate into pions, 

m,(Td) = 2m, (Td), (14) 

due to the s-channel pole that occurs at this point, see Eq.(ll). For our set of parameters, 
we find Td = 198MeV, while Tm = 215MeV. 

In Fig. 3, we plot the scattering lengths a^^°'^ as a function of the temperature. At 
T = 0, we obtain the numerical values 

a° = 0.179 

and (15) 
= -0.047, 

in units of the pion mass. These are close to the Weinberg values (a^)^ = 7m^/(327r/^) = 
0.16 and (a^)^ = -2m^/(327r/^) = 0.044, as one would expect [4,6], with the difference 
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being attributable to the composite structure of the mesons. Experimentally, one has 
[17] — 0.26 ± 0.05 and = —0.028 ± 0.012, which are compatible with those from a 
later reference [18], a° = 0.20 ± 0.01 and a? = -0.037 ± 0.004. By comparison, chiral 
perturbation theory for SU(2) [3] gives a° = 0.20, and = -0.042. It is seen that 
our results underestimate the ChPT and experimental results. This may be so since our 
calculation contains only the lowest order terms in 1/Nc, and does not include any higher 
order physical processes, such as pion rescattering [2]. We nevertheless have the advantage 
over ChPT in that we are able to calculate our model results over the entire temperature 
range. 

Regarding the temperature dependence of the a^, we see firstly that a° varies only 
slowly with temperature until T ^ 140MeV, and then displays a steep singularity at 
T = T(i, due to the fact that this amplitude contains T4, or physically, that it contains the 
s-channel pole due to a exchange. It diverges yet again at the Mott temperature Tm, since 
it also contains the amplitudes Ti, T3 and T5. On the other hand, is seen to vary slowly 
with temperature for temperatures less than T ~ 200MeV. Since this scattering length 
also contains the amplitudes T3 and T5, it exhibits a divergence at the Mott temperature. 

We compare our calculation with the zeroth order forms for (a^)^ given earlier in 
Eq.(15), by blindly extrapolating them to finite temperature with the replacement m^r — > 
mjr(T) and — > fi^iT), as calculated in the NJL model. A direct calculation of at 
finite temperature [19] indicates that /tt is monotonically decreasing with temperature, and 
diverges at the Mott temperature. Thus one sees that the extrapolated Weinberg estimate 
follows the calculation for a? rather closely. On the other hand, it is unable to reproduce 
the singular structure at for aP . One should also note that this naive extension of the 
Weinberg formula to finite temperature makes a definite prediction, i.e. that the ratio 
aP{T)/o?'{T) = —3.5 is a constant, independent of temperature. This appears to hold at 
the 10% level for temperatures up to T ~ 145 MeV, and is thus an excellent approximation 
at low temperatures. It breaks down of course in the vicinity of for a°. We comment that 
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similar trends are observed at finite densities, but at T = [5], where a phase transition 
also occurs. 

The results shown in Fig. 3 have been obtained in a particular model, the NJL model, 
which has proven rather reliable at the temperature T = 0. One may ask the question to 
what degree the results in the neighbourhood of the phase transition can be trusted, in 
particular in view of the fact that the NJL model only describes the chiral aspect of the 
phase transition but not deconfinement. The feature that the scattering lengths diverge 
at the Mott temperature may have a simple geometrical origin. Physically, when the 
pion reaches threshold, its constituents become unbound, and the pion radius becomes 
infinite. This is expressed in the NJL model via the relation of the charge radius to i.e. 

oc 1/ < >, where < > is the mean pion charge radius squared. Using the Weinberg 
extrapolation, which according to our calculation is a good (but not perfect) approximation 
of our calculated results, we have that the \a{T)\ increase like \a{T)\ oc< > (T), with 
a = or a^. This increase of in the neighbourhood of the chiral transition is thus 

really a consequence of the pion becoming unbound, i.e. it incorporates implicitly the 
signal of deconfinement. Therefore the increase of T Tm has a rather clearly 

understood physical origin, and may thus survive any model dependence. However, the 
detailed structure of singularities as shown on Fig. 3 may reflect some particular features 
of the NJL model. 

We conclude this paper by speculating as to whether there is a possibility of ex- 
perimentally verifying the predicted singularities of the tt — tt scattering lengths in the 
neighbourhood of the chiral (and deconfinement) phase transition. The signal observed 
in the Hanbury-Brown-Twiss (HBT) experiment of tt — tt correlations is related to (a) 
the space-time structure of the pion emitting system, and to (b) the pion-pion final state 
interaction. Any anomaly observed in the HBT results may be a hint as to the importance 
of (b), or indirectly the n — n scattering lengths. One may perhaps even try to compare 
events with different pion temperatures. 
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Appendix 

The temperature dependent functions /(p), K{p) and L{p) are listed here for the 
kinematics at threshold, p = ipo,p) = {y/s/2, 0). One finds 

r dk e f{E) - f{-E) 



X f dk k'^ 1 8E'^ 



{A2) 



+ 2pEmf{-E)^-^i 



and 

[dkk^^_ 1 12^^ 64^;^ ■ 



In these expressions, / is the Fermi function, f{E) = [exp(/?£?) — 1 



(A3) 
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Figure Captions. 

Fig.l: Box and cr-propagation diagrams for tt — tt scattering. 

Fig. 2: T matrix amplitudes in units of the pion mass, shown as a function of temper- 
ature. The Mott temperature Tm and a dissociation temperature Td are indicated by the 
dashed vertical lines. 

Fig. 3: Scattering amplitudes and a^, in units of the pion mass, shown as a function 
of temperature. The Mott temperature Tm and a dissociation temperature are indicated 
by the vertical dashed lines. Experimental values (see text) at T = are also indicated. 
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